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Abstract 

We compute the three-loop contribution to the N = 4 super symmetric Yang-Mills planar 
four-gluon amplitude using the recently-proposed Higgs IR regulator of Alday, Henn, Plefka, 
and Schuster. In particular, we test the proposed exponential ansatz for the four-gluon 
amplitude that is the analog of the BDS ansatz in dimensional regularization. By evaluating 
our results at a number of kinematic points, and also in several kinematic limits, we establish 
the validity of this ansatz at the three-loop level. 

We also examine the Regge limit of the planar four-gluon amplitude using several dif- 
ferent IR regulators: dimensional regularization, Higgs regularization, and a cutoff regular- 
ization. In the latter two schemes, it is shown that the leading logarithmic (LL) behavior 
of the amplitudes, and therefore the lowest-order approximation to the gluon Regge trajec- 
tory, can be correctly obtained from the ladder approximation of the sum of diagrams. In 
dimensional regularization, on the other hand, there is no single dominant set of diagrams 
in the LL approximation. We also compute the NLL and NNLL behavior of the L-loop 
ladder diagram using Higgs regularization. 
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1 Introduction 



Recent years have witnessed significant advances in M = 4 supersymmetric Yang-Mills theory 
(SYM) in four dimensions. The AdS/CFT conjecture relating M = 4 SYM theory to maximally 
supersymmetric string theory on AdSs x S 5 reveals various integrable structures. One aspect of 
recent progress is a greater understanding of the structure of on-shell scattering amplitudes, at 
both tree level and loop level. For example, based on an iterative structure found at two loops [1], 
Bern, Dixon and Smirnov (BDS) conjectured an all-loop form of the maximally-helicity-violating 
planar n-gluon amplitude [2], which is believed to be correct for n = 4 and 5, but requires 
modification by a function of cross-ratios for six or more gluons [3-7] . 

The simplicity of the BDS form for four- and five-gluon loop amplitudes arises from a hidden 
symmetry of the planar theory, viz., dual conformal invariance [8]. Dual conformal symmetry is 
also present in the theory at strong coupling, and can be seen in a string theory description of 
the scattering amplitudes, which are identified with Wilson loop expectation values in a T-dual 
AdS space [9]. The dual conformal symmetry of scattering amplitudes is understood as the usual 
conformal symmetry of the dual Wilson loops. In fact, this scattering amplitude/ Wilson loop 
relation extends to weak coupling as well [3, 10, 11] (for reviews see [12, 13]). 

The dual conformal symmetry is anomalous at loop level due to ultraviolet divergences asso- 
ciated with the cusps of the Wilson loops (which correspond to infrared (IR) divergences of the 
scattering amplitudes). Anomalous Ward identities can be derived for the Wilson loop expecta- 
tion values [14], whose solution is unique up to a function of conformal cross-ratios. The BDS 
ansatz satisfies the anomalous Ward identity, therefore it is exact for Wilson loops with four and 
five cusps, since there are no cross-ratios in these cases. Assuming the Wilson loop/scattering 
amplitude duality, this implies the validity of the BDS ansatz for four- and five-gluon amplitudes. 

Dual conformal symmetry extends to a dual super conformal symmetry [15], which is a sym- 
metry of all tree-level amplitudes [15-17]. (This dual superconformal symmetry can also be un- 
derstood from the string theory point of view by means of fermionic T-duality [18-20].) The dual 
superconformal symmetry combines with the conventional superconformal symmetry of M = 4 
SYM theory to form a Yangian symmetry [21]. The IR divergences of the loop amplitudes, how- 
ever, a priori destroy both the ordinary and dual superconformal symmetries (and therefore the 
Yangian symmetry); this breakdown is explicitly seen in dimensional regularization. Unlike the 
dual conformal symmetry, the breaking of the ordinary conformal symmetry is not (yet) under 
control (see however refs. [22-24] for progress in this direction), and so the role of the Yangian 
symmetry for loop amplitudes is unclear. 

In practice, it is desirable to use a regulator that preserves as many symmetries as possible. 
Recently, Alday, Henn, Plefka, and Schuster (AHPS) presented a regulator which, in contrast to 
dimensional regularization, leaves the dual conformal symmetry unbroken [25]. In this approach, 
the SYM theory is considered on the Coulomb branch where scalar vevs break the gauge symme- 
try, causing some of the gauge bosons to become massive through the Higgs mechanism. Planar 
gluon scattering amplitudes on this branch can be computed using scalar diagrams in which some 
of the internal and external states are massive, regulating the IR divergences of the scattering 
amplitudes (for earlier references see [9,26-28]). The diagrams remain dual conformal invariant, 
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however, provided that the dual conformal generators are taken to act on the masses as well as on 
the kinematical variables. The diagrams that can appear in the scattering amplitudes are highly 
constrained by the assumption of (extended) dual conformal symmetry. There is one point in 
moduli space for which all the lines along the periphery of the diagrams have mass m, while the 
external states and the lines in the interior of the diagram are all massless. This is believed to 
be sufficient to regulate all IR divergences of planar scattering amplitudes. The original SYM 
theory is then recovered by taking m small. 

Using this Higgs regulator, AHPS computed the M = 4 SYM four-gluon amplitude at one and 
two loops. They showed that the planar two-loop amplitude satisfies an iterative relation analo- 
gous to the one that holds in dimensional regularization [1], and suggested that exponentiation 
might extend to higher loops in an analog of the BDS ansatz for the four-gluon amplitude: 



logM(s,f) = -\-y{a) 



log 2 (4)+ log 2 (±) -g (a) log(-L) +log(-* 



Vm 2 , 

+ \l(a) [log 2 (0 + vr 2 ] + c{a) + 0(m 2 ) (1.1) 

where M(s,t) is the ratio of the all-orders planar amplitude to the tree-level amplitude, 7(a) is 
the cusp anomalous dimension (3.3), and Qo(a) and c(a) are analogs of functions appearing in the 
BDS ansatz in dimensional regularization (3.5). Let us emphasize that the nontrivial content of 
the BDS ansatz is the statement about the finite terms; the IR singular terms of the amplitude 
are expected to obey eq. (1.1) (or eq. (3.5)) on general field theory grounds. 4 

In this paper, we explore whether eq. (1.1) continues to hold at three loops and beyond. 
We compute the three-loop four-gluon amplitude using Higgs regularization, assuming that only 
integrals invariant under (extended) dual conformal symmetry contribute. Mellin-Barnes tech- 
niques are used to evaluate the integrals, some parts of which are computed numerically. We 
numerically evaluate the results at a number of kinematic points, and obtain explicit expressions 
in several kinematic limits (e.g., s — t and also the Regge limit s >> t). In every case, our results 
confirm the expected exponential ansatz (1.1) at the three-loop level. 

An important difference between Higgs regularization and dimensional regularization is that in 
the former, IR divergences take the form of logarithms of m 2 whereas in the latter, IR divergences 
appear as poles in e, where D = 4 — 2e. A consequence of this is that, provided eq. (1.1) is 
valid, the L-loop amplitude in Higgs regularization may be computed by simply exponentiating 
log M(s,t) without regard for the 0(m 2 ) terms since they continue to vanish as m — > even 
when multiplied by logarithms of m 2 . In contrast, the BDS ansatz in dimensional regularization 
(3.5) specifies log M(s,t) up to terms that vanish as e — > 0. In exponentiating log M(s,i), these 
neglected terms can combine with the IR poles to give rather complicated contributions to the 
IR-finite L-loop amplitude. 



4 In rcf. [29] the transition from amplitudes in dimensional regularization to amplitudes where (part of) the 
IR divergences are regulated by (small) masses was investigated. In this way, the log 2 (m 2 ) and log(m 2 ) terms in 
eq. (1.1) can be understood as arising from a different multiplicative renormalization factor relative to dimensional 
regularization. It is also conceivable that by adapting the formalism of ref. [29] to the present case one could show 
that the terms finite as m 2 — > follow from the corresponding formula in dimensional regularization. This would 
imply that the BDS ansatz can be stated in a scheme-independent way. 
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To put the matter the other way around, in order to test eq. (1.1) one need not compute 
any 0{m 2 ) terms of the Higgs-regulated L-loop amplitudes because they cannot make any con- 
tribution to the IR- finite part of log M(s,t), whereas to test the BDS ansatz in dimensional 
regularization, one must compute 0(e) and higher terms in the lower- loop amplitudes to obtain 
all the IR- finite contributions to log M(s,t). This is one of several significant advantages that 
Higgs-regulated amplitudes have over their dimensionally-regulated counterparts. 

A second major focus of this paper is the Regge behavior of planar Af = 4 SYM amplitudes 
using Higgs regularization. One motivation for this is that it presents a different way of exam- 
ining the iterative structure of the theory, in which some results can be obtained to all orders. 
Equation (1.1) may be rewritten as 

logM(s,f) = -l 7 ( a )log(^) -So(a) log (-^) - Q Q (a) log (^) + ^ 7 (a) + 5(a) + 0(m 2 ) 

(1.2) 

in which the log 2 (s/m 2 ) terms have cancelled 5 to leave a single \og{s/m 2 ). Consequently, M(s,t) 
exhibits exact Regge behavior 

/ a \ a(t) — l 

M(s,t) = M(t,s) = f3(t){ — ) (1-3) 
where the all-loop-orders Regge trajectory is 

a(t) - 1 = ~ 7 {a) log(^) - Go(a). (1.4) 

The lowest order term of the trajectory, — alog(t/m 2 ), gives rise to the leading log (LL) behavior 
of the L-loop amplitude 

in the Regge limit, whereas higher-order terms in the cusp anomalous dimension contribute to 
the NLL (next-to- leading- log), NNLL, etc. pieces of the amplitude (see eq. (4.10)). 

There is a subtle question of order of limits that appears when considering the Regge limit. 
In fact there are two ways this limit can be taken. The first possibility (a) consists in taking 
the limit m 2 <C s,t first, and then taking the limit s»i. This order is implicit in eq. (1.2). 
The second possibility (b) consists in taking the limit s ^> t,m 2 first, and then taking m 2 <C t. 
A priori, it is not clear that the result will not depend on the way the limit is taken. Explicit 
calculation shows that no such ambiguity arises for the integrals appearing at one and two loops. 
At three loops, as we will show, the individual integrals give different contributions in the two 
limits; the three- loop amplitude M^ 3 )(s,t), however, is unchanged. 

It seems that both ways of taking the limit can be justified, with slightly different inter- 
pretations. 6 On the one hand, the Regge (a) limit seems more appropriate in order to make 



5 A similar cancellation occurs in the Regge limit using dimensional regularization [10,30-32]. Subleading-color 
corrections to the Regge trajectory in dimensional regularization were also considered in ref. [32]. 

6 There is an analogous question in dimensional regularization of whether to take e -> or s > t first [33], but 
see the Erratum in v5 [34]. 
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contact with results in the massless theory. On the other hand, the Regge (b) limit is natural 
for scattering amplitudes with masses (see ref. [35] and sec. 4.4 below). 

We show in this paper that the leading log behavior (1.5) stems entirely from a single scalar 
diagram, the vertical ladder, in the Regge (b) limit of the Higgs-regularized loop expansion. 
(This contrasts with dimensional regularization, in which ladder diagrams do not dominate in 
the Regge limit; the leading log behavior of the L-loop amplitude receives contributions from 
most of the contributing diagrams. This might have significant implications for recent discussions 
of multi-Regge behavior [33,36-40].) We also compute the NLL and NNLL terms of the L-loop 
vertical ladder diagram. Other Higgs-regularized diagrams also contribute to the NLL, NNLL, 
etc. terms of the amplitude, although we do not yet have an all-orders characterization of which 
diagrams contribute to each order in the leading log expansion. 

In ref. [41], the LL approximation to the gluon Regge trajectory was computed using a ladder 
approximation and an alternative IR regulator, in which the external lines are massless while all 
internal lines of the diagrams are given a common mass (which cuts off the IR divergences) 7 . In 
the LL limit, this "cutoff regulator" yields results identical to those found in this paper. 8 The 
reason for this, as we will see, is that the LL approximation to the vertical ladder is insensitive to 
the masses of the propagators constituting the rungs of the ladder, and this is the only difference 
between the cutoff and Higgs regulators. We also compute the NLL contribution to the L-loop 
vertical ladder diagram using the cutoff regulator, which differs from the Higgs regulator by 
scheme-dependent constants. 

The Regge (b) limit of the scattering amplitudes can be understood by performing the Higgs 
regularization at a different point on the Coulomb branch of the theory. In this case, the scattered 
particles are massive, whereas some of the internal lines remain massless, eliminating the collinear 
but not the soft IR divergences. We suggest how the Regge behavior of the four-point amplitude 
can be understood from the cusp anomalous dimension of a Wilson line with a non-light-like 
cusp. 

The paper is organized as follows: In sec. 2, we give a short review of dual conformal symmetry. 
Exponentiation of the planar four-gluon amplitude in dimensional regularization and in Higgs 
regularization are discussed in section 3, and the three- loop amplitude in Higgs regularization is 
computed. The Regge limit of the planar four-gluon amplitude in several regularization schemes 
is examined in sec. 4, and the first three terms in the leading log expansion of the L-loop vertical 
ladder diagram are computed. Most of the technical details are relegated to three appendices. 

2 Short review of dual conformal symmetry 

In this section we give a short review of dual conformal symmetry, with a particular focus on 
four-gluon amplitudes. Hints for dual conformal symmetry first appeared as an observation that 

7 This procedure is closely related to an off-shell regulator. If all physical states have a common mass m arising 
from a conventional Higgs mechanism, then external lines with p 2 ~ are off-shell from this point of view. 

8 The Higgs regulator considered in this paper does not seem directly applicable to computing the subleading- 
color contributions of TV = 4 SYM amplitudes, which are not dual-conformal- invariant, and which involve non- 
planar diagrams. The cutoff regulator might be more suitable in this respect. 
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the loop integrals contributing to planar four-gluon scattering amplitudes in Af = 4 SYM theory 
have special properties when written in a dual coordinate space [8]. 

Let us recall that the full four-gluon amplitude can be decomposed in a trace basis, the 
coefficients of which are referred to as color-ordered amplitudes [42-44] . The color-ordered planar 
(i.e., large N) amplitudes may be written in a loop expansion 

oo 

A(p i ,e i )=Y i a L AW(p i ,e i ) (2.1) 

L=0 

in the 't Hooft parameter [2] 

^0( 4 ^- 7E ) € (2-2) 

where 7e is Euler's constant, and e = |(4 — D) with e < to regulate IR divergences. In Af = 4 
SYM theory, loop corrections to the four-gluon amplitude have the same helicity dependence as 
the tree amplitude, so we factor out the tree amplitude to express the amplitude as a function 
M(s,t) of the kinematic variables s = (pi + p 2 ) 2 and t = (p 2 + P3) 2 only 9 

00 

A(pi,ei) = AWfaejM&t), M(s,t) = 1 + ^ a L M^ L \s, t) . (2.3) 

L=l 

At one loop, we have 

M^(s,t) = -h 1 (s,t) (2.4) 
with the one-loop integral l\ given by 

L (s t) = (e^a 2 Y I (Pi+P2) 2 (P2+Ps) 2 , , 

{ ^' J in / 2 k 2 (k + pi) 2 (k + pi+ p 2 ) 2 (k - p 4 ) 2 ' 1 ' 

where the external states are on-shell: p 2 = 0. After a change of variables to a dual coordinate 
space [8,45], 

Pi = x t ~ x i+i , x 5 = x 1 , (2.6) 

one obtains 

hM = ( e ^y f^ 2 ff , (2.7) 

where now the on-shell conditions read x 2 i+1 = 0. Note that this change of variables can be 
most easily done in a graphical way (see fig. 1). From eq. (2.7) one can see that for D = 4 the 
integral would be invariant under conformal transformations in the dual coordinate space [8,45]; 
hence the name "dual conformal symmetry." Due to infrared divergences one cannot set D = 4, 
and hence the aforementioned symmetry is broken, which is why such integrals were later called 
"pseudoconformal." It was found that at least up to four loops all integrals appearing in the 



9 In this paper, we follow the — h ++ metric conventions of ref. [25], so that s is negative for positive CM 
energy. The amplitude M(s,t) will be real for s and t both positive. 
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Figure 1: One- loop scalar box diagram representing the on-shell integral (2.5) together with its dual 
diagram (2.7). The numerator factor (pi + P2) 2 {P2 + P3) 2 = ^13^24 * s n °t displayed. 



four-particle scattering amplitude have this property [10,46]. All this hinted at some deeper 
underlying structure. 

From the practical point of view, assuming that only pseudoconformal integrals contribute to 
an amplitude proved to be a useful guiding principle 10 (see e.g. refs. [5,47]). On the other hand, 
the inevitable breaking of the symmetry for D 7^ 4 was not under control. This changed when it 
was realized that the (finite part of the) logarithm of the amplitude satisfies certain anomalous 
Ward identities, which were initially derived for Wilson loops and are conjectured to hold for 
scattering amplitudes as well. Assuming the Ward identities hold for the scattering amplitudes, 
they explain the correctness of the BDS ansatz for four and five scattered particles [3, 14]. 

Initially dual conformal symmetry could be applied to the n-gluon amplitude M n for maximally- 
helicity-violating (MHV) amplitudes only, as can be seen from the fact that the variables in 
eq. (2.6) do not carry helicity. In ref. [15], it was shown how to incorporate the helicity infor- 
mation and to define dual (super) conformal symmetry for arbitrary amplitudes, both MHV and 
non-MHV. The predictions of ref. [15] about how dual conformal symmetry is realized at tree 
and loop level (through an anomalous Ward identity) have by now been checked to one-loop 
order [48-51]. The status of dual superconformal symmetry at loop level, which is related to the 
conventional superconformal symmetry, is under investigation [22-24]. 

The symmetries mentioned above can also be seen at strong coupling using the AdS/CFT 
correspondence. In a groundbreaking paper [9], a prescription for computing scattering ampli- 
tudes at strong coupling was given. There, a bosonic T-duality was used that maps the original 
AdS space to a dual AdS space. Dual conformal symmetry can then be identified with the 
isometries of the dual AdS space (up to the issue of regularization). It was later shown that 
the bosonic T-duality can be supplemented with a fermionic T-duality [18-20], which leads to 
the counterpart of the dual superconformal symmetry found in the field theory. The analysis of 

10 Note that there can be subtleties about which integrals should be called pseudoconformal (see e.g. ref. [5]), 
having to do with the peculiarities of dimensional regularization/reduction. 



7 



ref. [18] is valid to all orders in the gauge coupling constant; however, just as in the field theory, 
introducing a regulator may break the symmetry. It would be interesting to address this question 
in the string theory approach. A somewhat related question is how the helicity dependence of 
the scattering amplitudes is encoded in the string theory setup. (At strong coupling, it is argued 
to be an overall factor, which can be ignored, but certainly this is not the case at lower orders 
in the coupling constant.) 

In a recent paper [25] , a new regularization of planar amplitudes inspired by the string theory 
setup of refs. [9, 18] was advocated (see also refs. [26-28]). In the string theory, besides the 
usual stack of iV branes at z = (z being the radial AdS coordinate, with the AdS radius 
normalized to unity), M further branes are placed at distances z\ = l/m^. The scattering takes 
place on the stack of M branes. In the field theory, this corresponds to going to the Coulomb 
branch of M = 4 SYM. Specifically, one starts with a gauge group U(N + M) and breaks it 
to U(N) x U(M) ->■ U(N) x U(1) M . This leads to masses \m h - m i2 \ for fields with labels in 
the U(M), masses \rrii\ for fields with mixed N,M gauge labels, while the U(N) fields remain 
massless. We then consider the scattering of fields with indices in the U(M) part of the gauge 
group and take N 3> M. (In other words, we drop all diagrams containing loops of particles with 
gauge indices M.) Therefore, the internal labels will all be in the N part of the gauge group, 
while the particles running along the perimeter of all Feynman diagrams will have mixed (M , N) 
labels, and hence will have massive propagators. The latter make the integrals IR finite, and 
there is no need to use dimensional regularization. Of course, strictly speaking one is considering 
a different theory, but the original theory is approached in the small mass limit. 11 

The string theory setup suggests that the planar amplitudes defined in this way should have an 
exact, i.e. unbroken, dual conformal symmetry. This is possible because there are now additional 
terms in the dual conformal generators that act on the Higgs masses. These additional terms 
come from the isometries of dual AdS space. At one loop, the integral (2.7) is replaced by 

T( s f d 4 x a (^ 3 + (m 1 -m 3 ) 2 )(^ 4 + (m2-m4) 2 ) 

I 1 (s,t;m 1 ,m 2 ,m 3 ,m 4 ) = / . 2 2v 2 2v 2 , 2v 2 , j\ , (2.8) 

J in 2 (x 2 a + m 2 )(x 2 2a + m 2 )(x 2 a + m 2 )(x 2 a + m 2 ) 

now subject to the on-shell conditions x 2 ^ i+1 = — (rrii — m i+1 ) 2 , with the identification m 5 = m\. 
As anticipated, I\ is annihilated by the extended form of dual conformal transformations, 12 



K 11 h = , (2.9) 



where 



dxf 



(2.10) 



11 Unfortunately, since the requirement of finiteness imposes keeping N 3> M, one cannot reproduce the non- 
planar scattering amplitudes of the original theory. 

12 This can be seen most easily by thinking of the masses rrii as a fifth coordinate of the dual coordinates, 
defining xf 1 = (x^,mi), and considering conformal inversions in this five-dimensional space. For further details, 
see ref. [25]. 
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Figure 2: Mass assignment of a generic planar diagram in the Higgs setup. The external lines correspond 
to on-shell particles with pf = — (m,j — mj + i) 2 . The internal dashed lines correspond to massless particles. 
We call particles of mass \rrii — mj+i| 'light' and particles of mass rrii 'heavy' since the former become 
massless when we consider the equal mass case rrii = m. 



Note that the integral is finite (for rrii ^ 0), and the symmetry is exact (i.e. unbroken), hence 
there is no anomaly term on the r.h.s. of eq. (2.9). From eq. (2.9) one can deduce that the 
functional dependence of 1\ is [25] 

I 1 (s,t;m 1 ,m 2 ,m 3 ,m4) = f(u,v), (2.11) 

where 

ra x my, m 2 m A 

u— — and v = — . (2.12) 

s + (mi — m 3 ) z t + (m 2 — m^y 

Similar restrictions hold for a higher number of external legs. This form of dual conformal 
invariance in the field theory was checked for a particular four-scalar amplitude at one loop [25]. 
There it was also shown that assuming this symmetry at two loops leads to an iterative relation 
similar to that which holds in dimensional regularization. 

In fig. 2, we illustrate where the Higgs masses appear in a generic n-point integral in this 
setup. Thanks to dual conformal symmetry, we can set all masses equal without loss of generality, 
in which case all external particles are massless, the particles in the outer loop of a diagram are 
massive, and all particles on the inside are massless. 

One can in principle write down all dual conformal integrals at a given loop order and for a 
given number n of external legs. The amplitude would then be given by a linear combination of 
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these integrals [10] 

M n = 1 + ^V (x) c(Z)Z, (2.13) 
x 

where the sum runs over all dual conformal integrals 13 X, with L{T) the loop order of the integral, 
and c(X) a (rational) coefficient (a number in the case of MHV amplitudes). Further restrictions 
result from the requirement that the diagrams must arise from a scattering process; for example, 
at two loops, the one-loop integral (2.8) cannot appear squared. Moreover, there are restrictions 
on the number of propagators and numerator factors. Finally, integrals that would be formally 
dual conformal invariant but are divergent despite the introduction of the Higgs regulator should 
not appear [10]. 



3 Exponentiation of the four-gluon amplitude 

In this section, we recall the exponentiation of the four-point amplitude of M = 4 SYM theory, 
first using dimensional regularization to regulate the IR divergences, and then using the Higgs 
regularization described in sec. 2. 



3.1 Exponentiation in dimensional regularization 

On the basis of ref. [1], Bern, Dixon, and Smirnov conjectured [2] that the dimensionally- 
regularized all-loop orders amplitude (2.3) satisfies 

oo 

log M(s, t) = J2 a£ [f (e) ( e ) M(1) ( s > * I ^) + c(£) + 0(e)] (3.1) 
e=i 

where M^(s,t) is given by eq. (2.4), e = |(4 — D), and 

f { %) = \^ + \*tG? + <?f? (3.2) 

with the cusp 7^ and collinear Q^p anomalous dimensions [52] given by 14 

00 

7 (a) = a ^ (i) = 4a - 4C 2 a 2 + 22C 4 a 3 + • • • (3.3) 
e=i 

Qo(a) = ga^ ) = -C 3 a 2 +(4C5 + yC 2 C3)a 3 + --- (3.4) 

13 Here we mean dual conformal integrals in the sense of the Higgs regulator of ref. [25], not the off-shell regulator 
of ref. [10]. 

14 We use the notation of ref. [2]. Note that 7(a) = 2 r cusp (a), where r cusp (a) is also widely used in the literature. 
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The constants and f!p and the 0(e) terms in eq. (3.1) are not known a priori. The BDS 
ansatz (3.1) may be re-expressed as 



1=1 



logM(M) =J2 
where [2] 

OO 



7 



4(£ef 



2£e 



c (a 



V, 



2 +cw 



2\ 



7 (a) 



log 2 



120 



341 



17 



/• _ ^ 

216 9 3 



+c(a)+C(e) 
(3.5) 

••• (3.6) 



Overlapping soft and collinear IR divergences are responsible for the 1/e 2 pole in eq. (3.5). 

While the BDS ansatz (3.5) implies that the IR-finite part of the logarithm of the amplitude 
is simply expressible in terms of log(s//i 2 ) and log(t//i 2 ) and a set of constants 7^, Qq \ and 
the same is not true of the L-loop amplitudes M^ L \s,t) themselves. For example, eq. (3.5) 
implies [1] 



M^(s,t) = i[M«(M)] 2 - 



7 



(2) 



(2) 



+ 



Q(2) 
»0 



log 



/// 



+ log 



8e 2 2e 
t 



+ 



7 



(2) 



/// 



7 



(2) 



log 



/// 



+ log( 



(3.7) 



log 



/// 



log 



7T 



+ 7T7 



(2) 



+ c< 2 > + 0(e) . 



Because of interference between the positive and negative powers of e in [M^(s, t)] 2 , the 0(e~ v ) 
and O(e ) terms [53] in M^ 2 )(s,t) depend on more complicated functions (polylogarithms) of s 
and t, which are present in the 0(e) and 0(e 2 ) terms [1] of M^(s,t). In general M^ L \s,t) will 
receive contributions from the coefficients of positive powers of e in all lower-loop amplitudes. 



3.2 Exponentiation in Higgs regular izat ion 

The Higgs mechanism reviewed in section 2 can be used as a gauge- invariant regulator of the IR 
divergences in the planar massless theory, with the IR divergences appearing as log(m 2 ) terms 
in the amplitude, and any terms that vanish as m — > are dropped. One could ask whether, 
when regulated in this way, the four-point loop amplitude satisfies iterative relations similar to 
the BDS ansatz for the dimensionally-regulated amplitude. In ref. [25], it was suggested that the 
analog of eq. (3.5) is 



log M(s,t) 



~ g7(a) 



+ g7(a) 



log : 



log" 



(-) 

(D 4 



+ log 2 



+ c(a) + 0(m 2 ) 



(3.8) 
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where 7(a) is the cusp anomalous dimension (3.3), and Qo( a ) an d c(a) are the analogs of Qo(a) 
and c(a), but need not be the same functions since they are scheme-dependent [25]. Overlapping 
soft and collinear IR divergences 15 are responsible for the double logarithms in eq. (3.8). 

The L-loop amplitudes M^ L \s,t) are obtained by exponentiating eq. (3.8). In contrast to 
dimensional regularization, there is no interference between the IR-divergent log(m 2 ) terms and 
the 0(m 2 ) terms in eq. (3.8) since such terms vanish for m — > order by order in the coupling 
constant. Hence the amplitudes are simply expressed in terms of products of log(s/m 2 ) and 
log(t/m 2 ) and a set of constants 7^, Qq \ and c^. For comparison with later calculations, we 
explicitly write the predictions for the first few loop amplitudes, with u = m 2 / s and v = m 2 /t, 

M« = 
= 

+ 



+ 
+ 

where we have explicitly set 7^ = 4 and = = 0. 

The dual-conformal integrals that contribute through two loops are (see figs. 1 and 3) 

2 

M(s,t) = 1 - ^h(s,t,m 2 ) + °- [J 2 (s,t,m 2 ) + I 2 (t,s,m 2 )] + C(a 3 ) . (3.12) 

These integrals were computed in ref. [25], and the exponential ansatz (3.8) was verified to two- 
loop order. To determine the values of the constants in eqs. (3.9) and (3.10), it is sufficient to 
evaluate eq. (3.12) at s — t. Defining x = m 2 / s = m 2 /t, one finds [25] 

h(x) = 2\og 2 (x) -ti 2 + 0(x), (3.13) 
h{x) = log 4 (x)-^7r 2 log 2 (x)-4C 3 log(x) + ^7r 4 + C»(x). (3.14) 

15 In sec. 4.4, we will see that only soft divergences contribute in the Regge limit, giving a single logarithmic 
divergence. 
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-\og{v)\og{u) + --R 2 + 0{m 2 ), (3.9) 
X - \og 2 {v) log 2 ( M ) - Qtt 2 + \l {2 ^j log(u) log(u) 

Qf pog(«) + Iog( W )] + Qvr 4 + i7T 2 7 (2) + 5 (2) ) + 0{m 2 ) , (3.10) 

-\ \og\v) \og\u) + Qtt 2 + log» \og 2 {u) 

$j 2) [\og 2 {v) \og{u) + \og{v) \og 2 {u)] 
Qvr 4 + \ 2 ^ + + i 7 (3) ) log(«) log(u) 

Qtt 2 ^ + pog(«) + Iog(«)l 

(z^ 6 + iV 47(2) + l* 2 ^ + l n2l(3) + 5(3) ) + ° {m2) ' (3 - n) 



h ha hb 

Figure 3: Two- and three- loop four-point dual conformal diagrams. I2 and I 3a are ladder diagrams; 
is the tennis-court diagram. Numerator factors (including a loop-momentum-dependent factor for the 
tennis court) are omitted. See eqs. (A. 20) and (A. 23) for explicit formulae for the integrals. 




Figure 4: Three-loop four-point dual conformal diagrams, with numerator factors omitted. Both 
integrals require a factor of m 2 in the numerator in order to be dual-conformally-invariant. 

As discussed in ref. [25], these are consistent with eqs. (3.9) and (3.10) provided 

7 (a) =4a-4C 2 a 2 + 0(a 3 ), g (a) = -( 3 a 2 + C(a 3 ), 5(a) = — a 2 + £>(a 3 ). (3.15) 

A. ZAj 

The expression for 7(a) is consistent with eq. (3.3). 

We now test exponentiation (3.8) at the three-loop level, i.e., the prediction (3.11). Since 
an explicit calculation of higher-loop amplitudes using the Higgs regulator is not (yet) available, 
we will start from the assumption (2.13) that only dual-conformal integrals contribute. Equa- 
tion (2.13) requires two ingredients: the set of dual-conformal integrals and their coefficients. In 
order to identify the allowed set of dual conformal integrals it is helpful to use the dual notation 
and graphs introduced in section 2. (For more details and examples, see refs. [8,10,46,47,54].) 

There are four dual conformal integrals that can in principle appear [54] in a four-point 
three-loop amplitude: I 3a , hb, he and hd- The first two, depicted in fig. 3, are natural dual 
conformal analogs of the integrals appearing in the three-loop amplitude computed in dimensional 
regular izat ion [2]. The last two, depicted in fig. 4, are absent in dimensional regularization. We 
will assume that these are the only integrals required in the Higgs regularization, with the same 
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coefficients as in the dimensional regularization result. In general it is not valid to take a result 
computed in one regularization and transpose it to a different regularization; here this procedure 
can be justified a posteriori (as we discuss below) by imposing that the IR singular terms of the 
amplitude obey the relation eq. (3.8) as required on general field theory grounds. 16 
Given these assumptions we write 

M^(s,t) = ~ [l 3a (s,t,m 2 )+I 3a (t,s,m 2 ) + 2I 3b (s,t,m 2 ) + 2I 3b (t,s,m 2 )} . (3.16) 
o 

In appendix A, Mellin-Barnes (MB) representations for the integrals appearing in eq. (3.16) are 
derived. The small m 2 limit of these MB integrals is extracted using the same method as in 
ref. [25]. At the kinematic point s — t, we find 

ha{x) = f Q \og\x) + \n 2 \og\x)-\^oi^) (3.17) 
-64.93939402 log 2 (x) - 200.29103 log(x) - 196.597 + 0(x) , 

AO o 8 

h b (x) = — \og G (x)--n 2 \og 4 (x)--( 3 \og\x) (3.18) 

+37.8813132 \og 2 (x) + 113.11769 \og(x) + 90.0915 + 0(x) , 

where the decimal coefficients are approximations obtained by numerical integration of MB in- 
tegrals. We find that the three- loop amplitude (3.16) at s — t is consistent with the exponential 
ansatz (3.11) provided that 

7 (3) « 23.81111114+ 10~ 8 , ^ 3) « 2.68887+ 10" 5 , c (3) « -9.249 ± 10~ 3 . (3.19) 

together with eq. (3.15). We note that 7^) 22^4 is indeed the correct three-loop cusp anomalous 
dimension (3.3). 

Having obtained the coefficients from the s — t case, we are now in a position to test the full 
consistency of eq. (3.16) with eq. (3.11). In order to do this, we have evaluated the coefficients 
of the small m 2 expansion of the integrals I 3a and I 3b numerically for various values of the 
kinematical variables s, t. We have found agreement with eq. (3.11) within the numerical accuracy 
of the calculation. 

We assumed above that the coefficients of the integrals I 3c and I 3 d vanish in Higgs regular- 
ization, as they do in dimensional regularization. We have nevertheless evaluated these integrals 
in the small m 2 limit using MB methods, and found that 

I 3c (u,v) = 56.23 + 0(m 2 ) (3.20) 
I 3d (u,v) = -17.32 log(u) -62 + 0(m 2 ) (3.21) 

so that, if present, they could be simply accommodated in the BDS ansatz at three loops by 
redefining 17 £?q 3 ^ and in eq. (3.11). We will therefore not discuss them further in this paper. 

16 It would be desirable to determine the coefficients of these integrals using a unitarity-based method. Indeed, 
once a basis of integrals has been established, in our case using dual conformal symmetry, (generalized) unitarity 
cuts are a powerful tool to compute the coefficients of the integrals (see e.g. refs. [47,55-58]). 

17 Such a change could be detected when checking eq. (3.8) at higher orders in perturbation theory. 
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Equation (3.11) predicts the three-loop amplitude for arbitrary values of u and v. We would 
like to present a further test of eq. (3.11) in the limit s»i, i.e. This is Regge limit 

(a) discussed in the introduction. In order to find formulas for our integrals in Regge limit (a), 
we perform the small m 2 and subsequently the it <C ti limit in the MB integrals. In this way, 
we obtain an expression in terms of powers of log(w) and log(f), whose kinematic- independent 
coefficients are either numbers or (relatively simple) MB integrals. Where necessary we evaluate 
the latter numerically. 

Let us first collect the results for the one- and two-loop diagrams [25]. The one- loop box 
diagram gives 



lim lim h(u, v ) = log(w) 

«C» «,ti<l 



21og(u) + 0(v) 



+ 



IT 2 + 0(V) 



+ 0{u) . 



(3.22) 



For the two- loop horizontal ladder diagram (cf. fig. 3), one has 



lim lim I 2 (u,v) = log(w) 

ll<Cl> «,tl<l 



+ 



■\ \og\v) - 27T 2 l0g» - 4C 3 \0g{v) - -^TT 4 + 0(V) 

3 15 



(3.23) 
+ 0{u) , 



while for the two-loop vertical ladder 



lim lim I 2 (v,u) = log 2 (w) 

+ log(u) 

"1 



2\og 2 (v) + 0(v) 
-Uog 3 (v)-K 2 \og(v)-4( 3 + 0(v) 



(3.24) 



+ 



^og\v) + 27c 2 log 2 (v) + -7r i + 0(v) 



+ 0{u) . 



Substituting these results into eq. (3.12), we find agreement with eqs. (3.9) and (3.10) with the 
coefficients given in eq. (3.15). (The limit u«w does not affect the form of eqs. (3.9) and (3.10).) 

We now turn to the three-loop diagrams. For the three-loop horizontal ladder diagram (cf. 
fig. 3), we obtain 18 



lim lim I 3a (u,v) — \og(u) 

U<^.V li,t)<l 



^\oi{v) + j^log 3 ^) + gvrMog^) + 0(v) 



+ 

+ 0(u) , 



-Z-log» _ Z^log 4 ^) - ^Calog 3 ^) - ^7r 4 log 2 (t;) - 35.786 log(v) - 323.7 + 0(v) 
9U 9 o 4o 



(3.25) 



18 The coefficients of the it 4 terms were obtained numerically, and replaced by their probable rational equivalents. 
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while for the vertical three-loop ladder, we obtain 



lim lim I 3a (v,u) = log (u) 



90 + 



+ log 4 (w) 

+ log 2 («) 

+ log(«) 
2 



+ log>) 
44 



+ log 5 (w) 
4 



-— log(u) + 0O;) 



16 

^log 3 (t;)-— tt 2 log(u) + 0(u) 



+ 



2vr 2 log 2 ^) - 8(3 log(u) + ^tt 4 + 
2 74 

— log 5 (t;) + 8C3 log 2 (t;) - -7r 4 log(u) + 75.717 + 0(v) 
1 



(3.26) 



-^log» - -7r 2 log» - -( 3 \og 3 (v) - 111.5 log(u) + 141.2 + 0(v) 



+ 0{u) 



For the tennis court diagram in the orientation shown in fig. 3 we find, 



lim lim I 3b (u,v) = log (u) 



1 

90 



0(v) 



+ log 4 (w) 

+ iog 2 H 

+ log(w 

1 



-\\og\v)- 2 -^ + 0{v) 



+ \og 6 {u) 
AA 



+ log 5 (w 
16 



15 



log(u) + 0(v) 



1 6 

Q log 3 (t;) + — 7T 2 log(t;) + 0(t;) 



-- \og\v) - 4tt 2 log 2 (t;) - ^tt 4 + 0(v) 
o 45 

1 20 20 

- log 5 (w) + — 7T 2 log 3 (t;) - 8C3 log» + — 7T 4 log(u) - 24.8863 + 0(v) 
3 9 9 



(3.27) 



+ 



16 X3 

m \og 6 (v) + yCslog 3 ^) - -7r 4 log 2 (t;) + 111.499 \og(v) - 206.1 + 0(v) 



+ 0{u) , 



while for the tennis court in the transposed orientation, we obtain 



lim lim I^{v,u) = log 6 (w) 

U<^V U,t)<l 



180 



+ log>) 
+ log 2 (w) 
+ log(u 

1 



l -\og\v) + 1 -* 2 + 0{v) 



+ 0(v) 

+ log 3 (w) 
22 



+ log 5 H 



-— log^ + O^) 



-log 3 (t;)--7r 2 log(v) + 0(v) 



-\og\v) + -7t 2 \og 2 (v) + -7r i + 0(v) 



-^log 5 (v) - Vlog 3 (t;) - VlogM + 0(v) 



(3.28) 



+ 



log 6 (t;) + ^tt 2 log 4 (t;) - |Ca log 3 (v) + ^tt 4 log 2 (t;) - 24.888 log(u) + 280.8 + 0(v) 
9 3 15 



18 

+ 0(u). 

Plugging these into eq. (3.11), we see that the log 6 (w), log 5 (w), and log 4 (w) terms cancel, and that 
the log 3 (w), log 2 (-u), log(w), and log°(-u) terms of eq. (3.11) are reproduced using the coefficients 
given in eqs. (3.15) and (3.19). 
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In summary, we used the assumption of extended dual conformal symmetry to restrict the 
integrals allowed to appear in the three-loop four-point amplitude to I 3a , I 3 b, I 3c , and I 3 d- We 
then successfully checked the consistency of this assumption with the exponential ansatz (3.11). 
The first test consisted in computing the small m 2 expansion of the amplitude, and numerically 
comparing to eq. (3.11) for various values of the kinematical variables u and v. The second test 
consisted in comparing our result to eq. (3.11) in the limit u <^ v. 

4 Regge limit of the four-gluon amplitude 

In this section, we examine the four-point amplitude M(s,t) in the Regge limit s ^> t > 0. 
(Recall that M(s,t) is real in this kinematic region (see footnote 9). The "physical" Regge limit, 
— s ^> t > 0, differs from this by an imaginary contribution.) First we review the situation for 
dimensional regularization, and then we reexamine the Regge limit using Higgs regularization 
and cutoff regularization. The divergence structure of the Regge limit can also be understood 
via the scattering of massive particles. 

4.1 Regge limit in dimensional regularization 

We recall that in dimensional regularization the Regge form of M(s,t) is [10,30-32] 

/ e \ ce(t) — 1 

M(M) = /3(f) (f) (4.1) 

where the trajectory function is 19 

oo 

a(t)-l = £V 

e=i 

and the form of the residue /3(t) may be found in ref. [30,32]. 

The form of eq. (4.1) suggests that the IR-finite contribution to the L-loop amplitude (s, t) 
grows as log L (s) in the Regge limit. This conclusion, however, is not quite warranted because 
the neglected 0(e) and higher terms could in principle conspire with the pole terms to yield 
higher powers of log(s). Nevertheless, it was shown that this does not occur at least through 
three loops [32], and that the L-loop amplitude in fact grows as log L (s). It is not unreasonable 
to hope that this should hold to all orders. 

On the other hand, it is definitely not true that the individual L-loop diagrams contributing 
to M( L \s,t) are bounded by log L (s) in the Regge limit. As an explicit counterexample, consider 
the Regge limit of the contributions to the three-loop amplitude, keeping only terms of order 

19 The tree amplitude supplies another power of s/t so that, in the Regge limit, the color-ordered amplitude 
A(pi,ei) goes as (s/t) a(t K 



Me) 
Me 2 



ie 



+ 0(e) 



(4.2) 
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log 3 (s/t) or higher: 



ha(s,t) 
ha(t, S) 

hb(s,t) 
hb(t, s) 



0(C), 

t 



(4.3) 



-(—C 5 + —C 3 ) + O(e ) 
e \60 36 > 



2\ 3e 



12 



60 



13tt 2 
36 



0(6°) 



2\ 3^ 



6 



1 



\ 24^ J + e 



1 c^'-^c 3 



120 



72 



0(6°) 



0(£ 2 ), 

0(C 2 ) 



where we have denoted £ = log(s/t). Combining these contributions using eq. (3.16), we obtain 
the Regge limit of the three-loop amplitude 



M^(s,t) = 



3c 



1 

6^3 



7T 

24e 



+ 0(e c 



C 3 + 0(C 2 ) 



(4.4) 



but it is clear that the Regge behavior comes from no single diagram, and in fact a cancellation 
of higher powers of \og(s/t) occurs among the various contributing diagrams. This means that 
in dimensional regularization, it is very difficult to ascertain the Regge behavior of M^ L \s,t) by 
considering the Regge behavior of individual diagrams. 

In the next section, we will compare this to the Regge limit in Higgs regularization. As was 
noted in the introduction, there is a subtlety, namely that the Regge limit of individual diagrams 
can depend on the order in which the various limits are taken. As we will see in the next 
section, with one way of taking the limits, the leading Regge behavior of the L-loop amplitude 
is determined by a single contributing diagram, the vertical ladder diagram. 



4.2 Regge limit in Higgs regularization 

The conjectured exponentiation eq. (3.8) of the four-point amplitude in Higgs regularization can 
be rewritten as 



log M(s, t) = - ly(a) log (±) log ( JL) - So (a) 



7T 



+ '—l(a) + c(a) + 0{m 2 ) 
(4.5) 

As pointed out in sec. 3, in Higgs regularization we may simply exponentiate eq. (4.5) to obtain 
the four-point amplitude. Moreover, eq. (4.5) implies that the amplitude is Regge exact, 20 up to 
terms that vanish as m 2 — > 0: 



M(s,t)=f3(t) 



S \ <*(*)-! 



+ 0{m 2 ) 



(4.6) 



20 



Equation (4.6) could also be rewritten with s and t exchanged, due to the sot symmetry of eq. (4.5). 
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with 



a(t)-l = -77(a) logf-^ 



rrr 



P(t) = exp 



G (a) log + ^ 7 (a) + 5(a) 
V?7r / 8 



(4.7) 
(4.8) 



The — | 7 (a) log(t/m 2 ) piece of the trajectory a(t) agrees with the conjecture in eq. (2.18) of 
ref. [41], but in addition there is a contribution from the analog of the collinear anomalous 
dimension. 

Equation (4.5) implies that the leading-log (LL) behavior of the L-loop amplitude is given by 



M^(s,t) 



(4.9) 



In fact we can use eq. (4.5) to organize the L-loop amplitude in a leading-log expansion in 
log(s/m 2 ), with LL, NLL, and NNLL terms given by (recalling that u = m 2 / s and v = m 2 /t) 



M< L > 



'-Xogv) 1 



\og L {u) 



+ 



+ 



+ 



+ 



71 



+ 



7 



(2) 



2{L-l)\ 4(L-2)! 

C (2) 
^0 



(_l 0g ^-i + 



(L-2)! 



\ogv 



,L-2 



\og L -\u) 



(L-2)! 



(-logu) 



L-l 



' 8g (2) + ^ + 7r 2 7 (2) 27 (3) + 7r 2 7 (2) ^(2)^,2 



8(L-2)! 



'2^ 3) + vr 2 ^ + 7 ( 2 )^ 2) 



8(L-3)! 32(L-4)! 

^(2)x2 



logf 



,L-2 



2(L-3)! 4(L-4)! 
+ G(log L - 3 ( M )) 



L-4 



2(L-4)! 



log L - 2 (w 



(4.10) 



where we have explicitly set x ' = 4 and Qq = &• > = 0. In fact we have already tested the 
validity of this expansion through three loops in sec. 3 by evaluating the integrals that contribute 
to the amplitudes in a leading log expansion. 

As discussed in the introduction, however, there is a subtlety in how one evaluates the leading 
log expansion of the integrals contributing to the four-point amplitude. One approach, which we 
term (a), is to first evaluate all the integrals in the small m 2 limit, and subsequently to take the 
limit s t. This is what has been done for all integrals up to this point in the paper. In this 
section, we will explore another approach, which we call (b), in which we evaluate the s>i limit 
of the integrals first, for finite m 2 , and subsequently take the small m 2 limit of the coefficients of 
the leading log expansion. We will find that, although the intermediate details differ, the final 
result is still given (at least through three loops) by eq. (4.10). 
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For one and two loops, we have verified that both methods (a) and (b) yield the same results 
for the integrals namely, eqs. (3.22), (3.23), and (3.24). At three loops, however, the results begin 
to differ. The three- loop horizontal ladder diagram (cf. fig. 3), gives 21 



lim lim I 3a (u, v) = log(u) 

tt,D<Cl «<tl 



A log 5 ( w ) + j^ 2 log 3 ^) + gvr 4 \og(v) + 0(v) 



(4.H) 



+ 

+ 0(u), 



1 log 6 (u ) _ ^ 2 log4(u) _ log3(u) + 58^4 log2(u) _ 35 78Q hg{v) _ 323 ? + 



which in fact agrees with eq. (3.25). The vertical three-loop ladder, however, gives in the (b) 
limit 



lim lim h a (v , u) = log 3 (u) 



(4.12) 



+ log 2 (u) 

+ log(u 
34 



10 



-- log 4 (t;) - -Ti 2 log 2 (v) - 8C 3 log(u) + 0(v) 

^ log 5 (t;) + ^vr 2 log 3 (t;) + 8C 3 log 2 (t;) + ^tt 4 \og(v) + 75.717 + 0(v) 
Id 9 to 



+ 

+ 0{u). 



~ log» - |V log» - ^3 log 3 (t;) - ^tt 4 log» - 111.504 log( W ) - 363.4 + 0(v) 



In the other way (a) of taking the limit, the leading term goes as log 6 (w). For the tennis court 
diagram in the orientation shown in fig. 3 we obtain 



lim lim I^ b (u,v) = log 2 (u) 

u,v<^.l tt<Cf 



Uog\v) + ^ 2 \og 2 (v) + 0(v) 



(4.13) 



log(u 
"43 



q 44 7R 

-- log 5 (t;) - -7T 2 log 3 (t;) - 8C 3 log» - -tt 4 log(v) - 24.886 + 0(v) 
5 9 45 



Q9 16 14S 

^log 6 (v) + -7r 2 log» + -C 3 log 3 (t;) - — 7r 4 log 2 (t;) + 111.499 log(v) + 298.5 + 0(v) 
60 9 3 45 



+ 

+ 0(u), 

while for the tennis court in the transposed orientation, we obtain 

16 



lim lim hb(v , u) = log(-u) 



_ \og 5 (v) + -tt 2 log 3 (t;) + ^tt 4 log(^) + 0(v) 



(4.14) 



+ 

+ 0{u) . 



-A \ og \ v ) - \ og \ v ) - l& l g 3 (t;) - gvr 4 log 2 (t;) - 24.888 log(v) + 28.5 + 0(v) 



21 The coefficients of the it 4 terms were obtained numerically, and replaced by their probable rational equivalents. 
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Thus, the form of the leading log expansion of the three-loop integrals depends on the order in 
which the limits are taken. 

Despite the different expressions for the individual integrals, we find that the three-loop 
amplitude itself is the same in both Regge (a) and (b) limits. Plugging the expressions above 
into eq. (3.16), we find that the expression (3.11) is reproduced, using the coefficients in eqs. (3.15) 
and (3.19), just as when we evaluated the integrals in the other order of limits. In sec. 4.4, we 
give an heuristic argument which helps to explain why the Regge limit of the amplitude should 
be independent of the order in which the limits are taken. 

A significant difference between these two approaches is that while, in the Regge (a) limit, 
most of the diagrams contribute to the leading log term (4.9) of the amplitude (as is also the 
case in dimensional regularization), in the Regge (b) limit, only one diagram, the vertical ladder, 
contributes to the leading- log behavior of the amplitude in the Regge limit. Using 

L 



Ilo,(v,u) + (all other L— loop diagrams) 



(4.15) 



one can show that the leading log behavior of the L-loop amplitude (4.9) can be wholly ac- 
counted for by the vertical L-loop ladder diagram Ilo,(v,u), suggesting that all other diagrams 
are subdominant (i.e., grow no faster than log i_1 (s)) in the Regge (b) limit. Using the methods 
of ref. [35], we derive in appendix C the LL, NLL, and NNLL contributions of the vertical L-loop 
diagram, obtaining 

L L 

log u log V 



lim lim 

2 L 



+ 
+ 



(L-l)! 
2 L 



L-2 



LI 



u 



(L - 1) log L+1 v + (L + 2)C 2 log w v + (L - 1)C 3 log^ 2 v 



L-l 



,L-2 



(L-2)! 
L(L - 2) 



log' 



u 



10L 2 - 14L + 3 
180 



l 0g ^^ + (^ + l) 2 7r 2 l0g L^ 



v + 



IT* log 1 



+ (L - 2 y- 3 w -% 



360 

f- e>(io g L - 3 u) . 



v + ^-r^(L + 1.826)tt 2 C 3 log^ 3 ' 



(4.16) 



It is clear from comparing eqs. (4.10) and (4.16) that other diagrams begin to contribute at the 
NLL order. It is difficult, however, to use the methods of ref. [35] to estimate the leading loga- 
rithmic growth for the other diagrams due to the presence of numerator factors in the integrals. 
(It would be interesting to extend the methods of ref. [35] to these cases.) We suspect that only 
a small subset of the allowed L-loop diagrams contribute to the NLL terms of M^; it would be 
nice to characterize which. 



4.3 Regge limit in cutoff regularization 

Mandelstam [59] has given certain criteria for establishing whether Reggeization occurs, i.e. 
whether an elementary field in a Lagrangian field theory lies on a Regge trajectory or not. It 
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depends on a "counting procedure," which must be carried out separately for each field of the 
Lagrangian, where for renormalizable theories one simply does the counting at j = 0, 1/2, or 
1 to consider the issue. It has been shown [60-62] that the elementary fields of renormalizable 
Yang-Mills gauge theories lie on a Regge trajectory if the theory has a mass gap, where theories 
with a Higgs mechanism provide a case in point. These arguments are purely local in nature, in 
that the trajectory of a gluon passes through j = 1 at the mass of the particle, but no global 
information about the trajectory is provided in this construction. 

In addition to the "counting criteria," a certain factorization condition among the helicity 
matrices of the scattering process must be satisfied if the Reggeization is to take place. This 
requirement results from the solution of the unitarity-analyticity equations satisfied by the scat- 
tering amplitude analytically continued in angular momentum. As a result an integral equation 
must be satisfied, where the potential V is the inhomogeneous term in the integral equation. If V 
is taken to be the kinematical-singularity-free, partial- wave projection of the Born approximation 
helicity amplitude, the solution to the integral equation provides an analytic form for the Regge 
trajectory which is equivalent to the LL approximation to the trajectory function. 22 

Given this background, in ref. [41] the Regge limit of the ladder approximation for the color- 
ordered, tree-approximation-stripped amplitude was considered, with cutoff regularization (i.e., 
all the propagators in the loop integrals are given a mass m) for gluon-gluon scattering in M = 4 
SYM. It was shown that 

where 

a(t) = l-alog(-^)+0(a 2 ). (4.18) 

Hence eq. (4.18) agrees with eq. (4.8) to lowest order in the coupling. The reason for this is 
that, as we saw above, (1) the leading log contribution to the L-loop amplitude in the Regge 
(b) limit in Higgs regularization comes from the vertical ladder diagram alone, and (2) as we 
show in appendix C, the LL contribution of the vertical ladder is independent of the masses of 
the propagators of the rungs of the ladder, which is the only difference between the Higgs- and 
cutoff-regulated diagrams. 

In ref. [41] it was conjectured that with the cutoff regularization 

a(t) = l-i 7 (a)log(^), (4.19) 

which satisfies 

a(m 2 ) = 1 , (4.20) 

as required by general principles. Note that eq. (4.19) coincides with eq. (4.8) from the Higgs 
regularization up to a scheme-dependent constant, which begins at 0(a 2 ). 

While Higgs and cutoff regulators yield identical results at the leading log level, one could ask 
whether they differ at the subleading log level. Consider the diagrams that contribute to the two- 
loop amplitude (3.12). Whereas in Higgs-regularization the middle line is massless, in the cutoff 

22 See sec. 2 of ref. [63]. 
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scheme, all internal propagators have a uniform mass m. 23 To compute the two-loop diagram 
with uniformly massive propagators, more MB integrals are needed relative to the Higgs-regulated 
case. For example, a naive use of AMBRE [64] yields a nine-fold MB representation. Using the 
methods described in appendix A, we have written down a seven-fold MB representation, from 
which we extract the following results in the Regge (b) limit m<w<1: 



lim lim J 2; cutoff (u, v) 
lim lim J 2; cutoff (v, u) 

tt,D<Cl tl<tl 



= log(u) 

= iog 2 (w; 

+ log(w) 



- \og\v) + -7T 2 \og{v) + 4C 3 + 0(V) 

2\og 2 (v) + 0(v 



-Uog 3 (v)-^ 2 \og(v) + ^( 3 + 0(v) 



+ 0(log°(u)) , 

(4.21) 

+ O(log (u)) . 



Comparing this with eqs. (3.23) and (3.24), we see that the (3 coefficients in the log(tt) terms are 
changed. There are presumably also changes in the uncomputed log°(w) terms. 

We can also compute the NLL terms of the cutoff-regulated L-loop vertical ladder diagram. 
The calculation of appendix C is unchanged, except that the middle rungs of the (2 + 2e)- 
dimensional diagrams depicted in fig. 9 become massive. As a result, the integral B 2 in eq. (C.25), 
for example, is changed to 



e~ 2e7E ^cutoff =41og 2 M + e 
as a consequence of which we have 



10 , 3 4 2, 4 / 

ylog U+ -7T l0g«- -C 3 



+ 



(4.22) 



lim lim I La . cutoff (y, u) = — log u\og v 
2 L 



■ log L_1 u 



(L-l) 
+ 0(\og L - 2 u) 

which should be compared with eq. (4.16). 



1 -{L - 1) log L+1 u + (L + 2)( 2 log 1 " 1 v-^(L-l)C 3 log 



L-2 



(4.23) 



4.4 Scattering of massive particles and Regge behavior 

In this section, we present a different approach which helps to explain why we obtained the same 
three- loop amplitude from eq. (3.16) in both Regge (a) and (b) limits, despite the fact that the 
contributing integrals (3.25-3.28) and (4.11-4.14) differ so drastically from one another. 

Consider the four-point amplitude on the Coulomb branch of M = 4 SYM theory, as reviewed 
in sec. 2, with the scattered particles satisfying the on-shell conditions 

p 2 = -(m, - m i+1 f . (4.24) 

23 Note that, in contrast with dimensional or Higgs regularization, it is not clear whether other integrals might 
not also contribute in cutoff-regularization, as could happen when one goes off-shell. 
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Figure 5: (a) One-loop scattering of massive particles, (b) Sample higher-loop diagram. Fat lines and 
thin lines denote particles of mass M and m, respectively. Dashed thin lines denote massless particles. 

Recall that the assumption of dual conformal symmetry implies that the amplitude is a function 
of two variables u and v, defined in eq. (2.12). In previous sections of this paper, we considered 
the equal mass case, rrii = m, so that u = m 2 / s and v = m 2 /t, and the external states are 
massless p 2 = 0. The Regge limit s > t > m 2 corresponds to u <C v <C 1. 

In this section, we consider instead the two mass case, mi = m 3 = m and 1712 = m± = M, 
which implies p 2 = — (M — m) 2 , and (cf. eq. (2.12)) 



u 



m 

— and 

s 



V = 



M 2 

T 



(4.25) 



Dual conformal symmetry implies that the Regge limit n<t)<l can be attained by choosing 
m 2 <C M 2 s,t, so that the on-shell condition becomes pf ~ —M 2 . This corresponds to the 
scattering of particles of mass M by the exchange of particles of much lighter mass m. See 
fig. 5(a) for the one-loop contribution to this scattering. At higher loops, massless particles will 
also be exchanged in the interior of the diagram; see fig. 5(b) for a sample higher- loop diagram. 

As was discussed in the introduction and in section 4 it is important to specify in which order 
the Regge limit is taken. In this section, we will first take m 2 small for fixed M 2 , s, and t. This 
corresponds to the first part of Regge limit (b), namely, u — > 0, i.e., u< « and m < 1. Later 
in this section, we will take M 2 s, t, which corresponds to the second part (v <^ 1). The 
amplitudes discussed above have soft IR divergences (as m — > 0), but no collinear divergences, 
so the double logarithms log 2 (m 2 ) characteristic of overlapping divergences are absent, leaving 
only single logarithms log(m 2 ). At one loop one can show that 



mW(u,v) — > --(io gM )r cusp ( 

u— >0 Z 



^Lioop + ^V) 



(4.26) 



where T cusp (a,9) is the cusp anomalous dimension [65-68] of a Wilson line with a non- light- like 
cusp, and with 9 being the (Minkowskian) cusp angle defined by cosh^ = p 2 ■ p 3 /M 2 . It is 
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plausible that the leading behavior for all loops is 

M(u,v) — > exp{-h\ogu)T cusp (a,9) + O(u )\ . (4.27) 

u — ^0 I Z I 

Next, we take the further limit M 2 <C t, i.e., ti « 1. In this limit, the cusp anomalous 
dimension goes to [69] 

r cusp (M) — > (logt;)r cusp (a) + 0(t; ), (4.28) 
d— >o 

where 6 1 — logv and r cusp (a) = r y(a)/2. Accepting the conjectured eq. (4.27), one obtains in 
the Regge limit 

lim lim M(u,v) = exp { -i(logw)(logt;)r cusp (a) + •••!. (4.29) 

Since we have taken the limits in the order tt< 1 with v fixed, followed by v <C 1, this is Regge 
limit (b) as defined in the introduction. Note that we arrive in this way at the same (leading log) 
result as eq. (4.5), which was obtained in Regge limit (a), i.e., taking m, t)< 1 first, followed by 
u <Cw. This helps to explain why we obtained the same three-loop amplitude from eq. (3.16) in 
the two Regge limits, despite the fact that the contributing integrals had different forms in these 
two limits. 

It is also noteworthy that eq. (4.29) directly exhibits the single log behavior expected of the 
Regge limit, without the cancellation of double logs that occurs in taking the Regge limit (4.5) 
of the exponential ansatz (3.8). 

5 Discussion 

In this paper we have tested various conjectures and issues related to the Higgs regulator scheme 
proposed in ref. [25] for J\f — 4 SYM theory. The assumption of exact dual conformal symmetry 
enabled us to compute the four-gluon scattering amplitude in the planar limit to three-loop 
accuracy. Our results are consistent, in various limits, with an analog of the BDS conjecture 
proposed for the Higgs-regulated four-gluon amplitude, with the same value (to three loops) of 
the cusp anomalous dimension as that found in dimensional regularization. This is as expected 
since the cusp anomalous dimension is IR-scheme-independent. 

The Regge limit of the four-point function was considered in various IR regulator schemes. 
It was shown that the leading log behavior of the Regge trajectory is determined by the sum of 
vertical ladder graphs for the Higgs regulator and for a cut-off regulator. 24 This contrasts with 
dimensional regularization, where no single set of diagrams dominates at any order of pertur- 
bation theory. Further one may associate the perturbative expansion for the cusp anomalous 
dimension with the NLL, NNLL, • • • approximations to the gluon Regge trajectory. A particular 
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Here we are referring to the Regge (b) limit, as defined in the paper. 
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dual conformal mass assignment allows one to obtain directly the single logarithmic behavior of 
the Regge trajectory, without the necessity of cancellation of log 2 terms. 

We also found in the course of this work that Higgs-regulated amplitudes lead to more efficient 
integral representations (i.e., fewer MB parameters) than more generic (e.g., cutoff- regulated) 
schemes. 

There are several issues deserving further attention. It would be useful to extend the methods 
of ref. [35] to determine the Regge limits of diagrams involving non-trivial (i.e., loop- momentum- 
dependent) numerator factors. It would also be interesting to determine whether there is a simple 
characterization of the subset of L-loop diagrams that contribute to the NLL, NNLL, • • • terms 
in the L-loop amplitude in the Regge limit. 

Overall, we have seen that the Higgs regulator for planar Af = 4 SYM amplitudes has a 
number of advantages over other regulators. 
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A Mellin- Barnes representations of the integrals 

In this section we derive Mellin-Barnes (MB) representations for the integrals considered in the 
paper. In doing so we follow the loop-by-loop approach advocated in refs. [72, 73]. In this 
approach, one successively derives MB representations for one-loop subintegrals. Remarkably, 
the Higgs masses can be incorporated rather naturally into this procedure. 

Tools that we have found useful are the Mathematica packages MB [70] and AMBRE [64]. 
The latter can also be used to derive MB representations. Of course, one can sometimes find 
MB representations involving fewer MB integrals by going through the derivation by hand. Since 
the number of integrals we need is rather small we have followed the latter strategy. We present 
details of the derivation below, since the essential steps when deriving MB representations for 
higher-loop or higher-point integrals are the same. 

It is interesting to note that using a uniform cut-off (i.e., giving a mass to all propagators of 
the integral) makes the resulting MB representations considerably more intricate, and one needs 
more MB parameters. 

A.l MB representation for the 3- loop ladder diagram 

Let us consider the three-loop ladder depicted in figure 6a, 

f d 1 ^ cc df^ cc q *tjC 'y , ^ 

ha = S t / , (P25,m Pl5,m P35,m P56 P36,m Pl6,m P67 P37,m Pn,m Pi7,m) ■ (A.l) 
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Figure 6: Diagrams of two dual conformal integrals at three loops. Thick lines correspond to 
massive propagators, thin lines to massless ones. The Xi are dual coordinates. Figure (a) shows 
the three- loop ladder diagram, with a factor of s 3 t = {x\^fx\ A removed. Figure (b) shows the 
tennis court diagram, with a factor of st 2 removed and a numerator x\ 7 + m 2 not shown in the 
diagram. 



Here Pij, m = x% + m<2 and Xi is the dual notation introduced in sec. 2. In the massless case, it 
is possible to derive a convenient MB representation by successively doing the loop integrations 
[72,73]. As we will see presently, the same strategy also works well in the massive case. We begin 
with the x 5 subintegral 

lia* = [ ~ ~T (Pl5,m P25,m P?,<o,m Pf>&) 1 (A-2) 



where 



ha = S 3 



7 . 2 \2 ^ a (^36: m Pl6,m P&7 P?>7,rn Pv7,m P&1 >) • (A-3) 



Introducing a parameters (see for example ref. [73]), we obtain 

/« = r mm (A . 4) 

Jo [aia 3 s + a e (aiP 16im + a 2 P26, m + a 3 P 36im ) + (a x + a 2 + a 3 ) 2 m 2 ] 



Note that the range of the sum in the delta function can be chosen arbitrarily. Here, it is 
convenient to choose X^ =1 in order to simplify the m 2 term. Performing the a§ integration, we 
obtain 

j(i) = f 1 d<*i <HE?=i a i - x ) / A 5 x 

J (^P^m + a 2 P 2 6,m + a3P36,m)(aia 3 s + m 2 ) ' 



27 



To perform the remaining integrations over the a parameters, we introduce three MB parameters. 
Using eqs. (A. 39) and (A. 40), we obtain 

i£ = |^r(^ 1 )r(-^)r(-^)r(i + ^ 1 + ^)r(i + ^)(m 2 )- 1 ^ 

3 

x / da^(^a s - 1) (a 1 a3^ 2 (a2-P2 6 , m )^ 1 " 1 " 3 (« 1 P 16 , m ) 2l («3i 5 3 6im ) 23 • (A.6) 
J i=i 

The parameter integrals are now easily done with the help of (A. 43), leading to 

4'J = /j^/^^sK^T (A.7) 

where 

f il \z 1 ,z 2 ,z 3 ) = T(-z 1 )T(-z 2 )T(-z 3 )T(l + z 1 + z 3 )T(l + z 2 ) 

r(-^ - * 3 )r(i + z 2 + ^)r(i + z 2 + z 3 ) 

r(2(i + %)) ' 1 ' ] 

Plugging (A.7) into (A. 3), we obtain 

ha = S 3 t J ^f^(z 1 ,Z 2 ,Z 3 )^(m')- 1 -^ I ^-4 2 J (P^P^mP^m)- 1 , (A.9) 

where if} is the integral over x§: 

4? = /^(P 26 , m )' 1 --- 23 (Pi 6 , m )- 1+ -(^ m )- 1+ -(P 67 )- 1 

don 5(Y, on - 1) «r 21 «3 23 «2 1+23 



Jo 



[a>ia 3 s + a 7 (aiPi 7 , m + a 2 P 2 7,m + "3^37,™) + («i + «2 + a 3 ) 2 m 2 ]' 



1 

x r(i-^)r(i-^)r(i + ^ + ^) ( ' 10) 

This is a generalization of the integral considered before with more general powers of the propa- 
gators. The calculation is completely analogous, and we obtain 



,-(2) 



I J2rif f(2){Zl > 2 > 3] Vw) ^ (m2r1 "" 2 {Pl7 > m)V1 {P ^ mr ^^~ l ~ Vl ~ m ' ( A - U ) 



with 



, {2)( , T(- yi )T(-y 2 )T(-y 3 )T(l + y x + y 3 )r(l + y 2 ) 

^ 1 > 2 ' 3 ' 2/1 ' 2 ^ " r(i-^)r(i-^)r(i + ^i + ^) 

r(i - + yi + y 2 )r(^i + * 3 - 2/1 - 2/3)1X1 - z 3 + y 2 + y 3 ) 

r(2(l + y 2 )) - [ ■ 1 
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At this point we note that, for a ladder diagram with more rungs, the next subintegral would be 
of the same type as the preceding one (up to a change of labels of MB parameters), see figure 7. 
In the next section, we will use this iterative structure to write down a MB parametrization for 
the L-rung ladder with L arbitrary. 

Coming back to the three-loop ladder, we are left with the final integration over x?, namely 



ha = S 3 



t J (A.13) 



where 
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4? = / ^ (P27,m)- 1 - yi - y3 (Pl7, m )- 1+yi (P 37 , m )- 1+ ^ (P^rn)- 1 ■ (A.14) 

Introducing parameter integrals we have 



(i - 2/3) ./„ w 



jyi+1/3 p-V3 

' ' r(i- yi )r(i + yi + y 3 )r(i-y 3 ) k + fafkt + (ft + J 2 + (3 3 + ft) 2 ™ 2 ] 2 ' 

(A.15) 

Obviously, here it is convenient to choose = S«=r Introducing two more MB parameters, we 
find 

4? = / / (3) (2/i,2,3; ^,5) S - (m 2 )— - , (A.16) 



where 

/ l3j (1/1,2,3; 24,5) = 



( 3 ), .„ , T(-z 4 )T(-z 5 )T(2 + z, + z 5 ) 



r(i- yi )r(i-y 3 )r(i + yi + y 3 ) 



r(l + z 5 )r(l +y! + y 3 + z 5 )T{l - yi + z 4 )T{l -y 3 + z 4 ) 
T(2(2 + Zi + z 5 )) ' 

Putting everything together we arrive at the final result (relabelling z 4)5 — > z 7t8 and 1/1,2,3 — > -24,5,6) 

T f 6^1,2,3,4,5,6,7,8 n ,( 2 )/ \ ,(3)/ \ ( ^ \ *° ^ ( ^ \ * 

ha = J 7^— y P J (^l,2, 3 ) P j (^l,2, 3 ; ^4,5,6) h ^4,5,65 Z 7 , 8 ) [~J J ^ — J 

(A.18) 

In eq. (A.18) the integration contours are chosen parallel to the imaginary axis, with the real 
parts of the integration variables defined such that all arguments of the Y functions have positive 
real part. An allowed choice is 

1 1 1 3 1 15 9 13 , . 

Zl ~ -~r Z2 ~ -~V Z '~ -~V Z '~ -~V Zh ~ -~W Z% ~ -~32' Z ' J ~ -~32' Zs ~ -"32- (A ' 19) 
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Figure 7: Diagrammatic representation of the iteration used to derive the L-loop MB represen- 
tation for the ladder integral. Fat lines correspond to massive propagators, thin lines to massless 
ones. The numbers indicate the powers of the propagators. One can see that after removing a 
rung of the ladder (by introducing three Mellin-Barnes integrals) we obtain an integral of the 
same type, and the procedure can be iterated. 

A. 2 MB representation for the L-loop ladder diagram 

As was noted in the previous section, one can straightforwardly write down an iterated MB 
formula for the L-loop ladder diagram. The decisive step is illustrated in figure 7. Let us define 



/L+4 , , 4 s L+3 
II(~^<f) (,P25,mPl5,mP35,m) * J\ {Pj,j+\Pl ,j+l,mPz,j+\,m) 1 (P<l,L+4,m) 



j=5 \ / j=5 

The MB formula we obtain is (3L — l)-fold and is given by (for L > 1) 



(A.20) 



Iia — 



L-l 



n l (2^3 ) (2^)2 ^ 8 



m 

T 



L-2 



For completeness, we also give a formula for L — 1, 

^1,2 u -i-,^-i-, 2 r(-^)r 2 (i + ^)r(-^ 2 )r 2 (i + Z2)V{2 + Zl + Z2 ) 



(A.21) 



ha = h 



(2vri) 2 ~ ' T(A + 2 Zl + 2z 2 ) 

which is equivalent but more symmetric than the one given in ref. [25]. 



(A.22) 



A. 3 MB representation for the tennis court diagram 

The procedure above can be repeated for the tennis court diagram (fig. 6b), which is given by 

d A x 5 d 4 x 6 d A x 7 



'36 



st z 



(«7T 2 ) ; 



Pn,m (P25,m Pl5,m L57 P56 Pl6,m PiQ,m P&7 Pl7 > P?>7,m Pi7,m) 1 • (A. 23) 
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A new feature of the tennis court is the presence of a numerator required by dual conformal 
symmetry. It is possible to treat this numerator as a propagator with negative power. However, 
this would require an analytical continuation, for example in the dimension of the integral, 
D — > 4 — 2e. Since our integral is finite in four dimensions, we find it preferable to perform the 
calculation in such a way that all steps can be done in four dimensions. This can be achieved by 
carrying out the two-loop ladder subintegral first (integrations over x 5 and x 6 ). The numerator 
will then combine in a natural way with a propagator obtained from this subintegral, allowing 
us to use a simple formula for the final integration over x 7 . 

As was explained above, we start by computing the subintegral 1$ defined by 



'3b 

j4„ 



where 



^sb ~ f ~~ir (-f25,m Pi5,m Phi P*>&) 1 (A. 24) 

J 27T Z 

hb = St 2 J Jj^Tj2 7 ^17,m lib (Pl6,rn P46,rn Pq7 P27,m P37,m P<l7,m) 1 ■ (A. 25) 



It can be written as 



_ f 00 da i 5{Y J ai + a 2 - I) 

[{PvT,m<X\ + P27,m0il)0L7 + (Pie^l + ^6,m«2)«6 + ^67«6«7 + ™ 2 \ 



(A.26) 



The integrals over a 6 and a 7 are done using eq. (A. 42). Introducing two further MB integrals to 
factorize both b and c in eq. (A. 42), we arrive at 

4 X) = / ^ (™ 2 r 1-21 4 1} (Pie, m r- z3 (p«, m r- n {p^r (p m )- i ~ n (A.27) 



with 

= r 2 (i + z 1 )r(-z 2 )r(-z 3 )r(-z 1 + z 2 )r(-z 1 + z 3 ) 

^ r(l +Z!+Z 2 - Z 3 )T(1 + Z1-Z2 + z 3 ) 

r(2(i + ^)) ' { } 

Next, we want to carry out the x 6 integration in eq. (A. 23). We define 

d 4 x 7 dz lj2 , 

in 2 (2m) 

(A.29) 

Collecting all propagators involving x 6 in eqs. (A. 23) and (A.27) we have 

4b = f^T (Psr)- 2 - Zl (Pi6, m )- 1+Zl - Z3 (P 26 , m ) 23 (P46,™)- 1 (A.30) 
J in 

1 f da i 8(^2 124 a i -l)a^ Zl+Z3 o^ 1 ~ Z3 aj +Zl 



2 r d X 7 dZi )2 fi „(1) , 2\-l-zi (p \l+z 2 (p \-l+zi-z 2 (p \-Wp \-l f 2 



T(2 + Zi)r(l - Z 1 + Z 3 )T(~Z 3 ) J [ta 2 a A + (Pi7 >m Q!i + P 2 7,mtt2 + P47,mOi4)a 7 + 



212 



; r / dCHiSiy OLi — 1) «i 

'-*) Jo ' ^ ' 



^( Zl ) I 7 w-V^ , N ---1+13^-1-23 v 



r(i - ^ + z 3 )r(- 

x (ta 2 a A + m 2 ) 21 (Pi 7 , m ai + Pvr, m a 2 + i^mt^) -2- * 1 • (A.31) 
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This is almost identical to an integral considered earlier for the ladder diagrams. Introducing 
MB parameters £4,5,6 we find 

4 2) = / fS ] (™ 2 r~ Z4 (p^r (P27, m r (p 47 , j- 2 -™- 26 , (A.32) 



with 



(2) T(-z 4 )T(- Zl + z 4 )T(-z 5 )T(-z 6 )T(2 + z x + z 5 + z 6 ) , 
36 T(2 + ^)r(l - Zl + z 3 )T(-z 3 )T(2(- Zl + z 4 )) x 

x T(l - Zl + z 3 + z 5 )T(-z 3 + z 4 + z 6 )T(-l - Zl + z 4 -z 5 - z 6 ) . (A.33) 



Finally, we carry out the remaining integration over x 7 , 

/£> = / ^ (P 17 , m ) 1+Z ^> (P27, m )- 1+Zl - Z2+Z6 (P37, m )- 1 (P47,m)- 3 - Zl - Z5 - Ze ■ (AM) 
J ITf 

Up to a relabelling, this is the integral ifj of eqn (A. 14). Hence we have 

7 » = j r ^-jt ) ^t z Hmr 2 -^ z \ (A.35) 



with 



,(3) r(-z 7 )T(-z s )T(2 + z 7 + z 8 ) 

T(l - yi )T(l - y 3 )T(l + Vl + y 3 )T(2(2 + z 7 + z 8 )) 

x T(l + z 7 )T{l + Vl +y 3 + z 7 )T{l - Vl + z 8 )r(l -y 3 + z 8 ) . (A.36) 



Here we have swapped s and t with respect to (A. 16) and y\ = —2—z\—z^—z & and y 3 = zi— Z2+-26- 
Collecting all factors we obtain our final expression for the tennis court diagram 

^1,2,3,4,5,6,7,8 (m 2 \ 1 Zr fm 2 \ 24 Zfi (x) ( 2 ) (3) 



An allowed set of real parts for the integration variables is given by 

11 1 7 1 9 3 13 1 

zi - - I g, ^2 - - g, ^3 - - jg, ^4 - -g, £5 - -g, ^6 - -33, ^7 - - 33, ^8 - ■ (A-38) 

A. 4 Auxiliary formulae 

Here we collect some auxiliary formulae that were used in the derivations above. It is understood 
that the formulae below should be used within their region of validity, i.e. the real parts of the 
arguments of the T functions should be chosen positively. 
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(a + b)-" = ^Y{-z 2 )V{\ + z 2 )a^b-^\ (A.39) 

1 C Ay Ay 

(a + b + c)- x = ^J a ^T{-z 1 )T{-z z )T{\ + z l + z z )a^b^c-^- z \ 

(A.40) 

dxx z ^a + bxr = a l+z ^b- 1 ^ F(1 + Zl)V }~ 1 ~ Zl ~ Z2) , (A.41) 

T(-z 2 ) 

l" (a + bx T^ +ixyY - (2«)- 1 /^M^c)-r^,)r=(i + ,), (A.42) 

/'n^d-fa,) = ^IM. (A,3) 



^0 



B Numerical evaluation using sector decomposition 

Here we present an alternative, semi- numerical approach for evaluating the small m 2 expansion 
of integrals that is independent of the various MB representations derived above. We use the 
availability of powerful numerical algorithms for the calculation of the e-expansion of parameter 
integrals [71,74]. The expansion is done analytically, and only the coefficients of the poles are 
evaluated numerically. Therefore, this method avoids the problematic large logarithms that 
would appear in a direct numerical integration. 

In order to use these algorithms, we need to reformulate our problem as a calculation of the 
e-expansion of some parameter integral. Consider a generic four- dimensional L-loop diagram 
containing a propagators, where propagators with indices i G M. have mass m 2 , and all other 
propagators are massless. 25 Then, the Feynman representation reads [73] 

2L—a 



I = T(a — 2L) J dai 5(1 - <*i) U a 



-2L-2 



V{s,t) + Um 2 ^2ai 



(B.l) 



i&M 

where U and V(s,t) are polynomials in the a« (see ref. [73] for more details). 

We introduce one Mellin-Barnes parameter in order to separate off the mass dependence, 

r dz 

1 = J —T(-z)T(a-2L + z)m 2z f(s,t,z), (B.2) 

f(s, t,z) = j da t 5(1 - «i) (E a ^ U a - 2L ~ 2+Z V(s, t) 2L - a ~ z , (B.3) 
J i ieM 

with 2L — a< Re(z) < 0. Note that f(s, t, z) is very similar to a dimensionally-regulated massless 
Feynman diagram of the same topology. 26 



25 Modifications are necessary for the case of integrals with numerator factors. 

26 A. Zhiboedov has independently written down a similar formula at one loop which exhibits the relation 
between massive and massless Feynman integrals (private communication). 
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Since we want to determine the behavior of / as m 2 — > it is convenient to shift the integration 
contour in z to positive values of Re(z). In doing so, one picks up a residue from the pole at 
z = 0, and it is clear that the logarithms in m 2 will come from (minus) this residue. 



-Residue [T(-z)T(a -2L + z)m 2z f(s, t, z)] 



+ 0(m 2 



(B.4) 



z=0 



The value of the latter can be obtained from the expansion of T(—z)T(a — 2L + z)m 2z and 
f(s,t,z) about z = 0. Therefore, in order to compute I up to and including order log°(m 2 ), we 
need to evaluate f(s,t,z) to order z°. Note that we do not need terms that vanish as z — > 0. 
The coefficients of this expansion can be computed numerically using ref. [71]. 

We have tested this method for the two-loop ladder and found agreement with the result of 
ref. [25]. The method is also applicable to three- loop integrals but requires more computer time 
and memory. 



C Regge limit of the L-loop ladder diagram 

In the two- and three-loop amplitudes considered in this paper, we see that the leading log 
behavior in the Regge (b) limit (see discussion in the introduction and in sec. 4) arises from the 
vertical ladder diagram. In this appendix, we will compute the first three leading terms in the 
Regge limit (b) of the ladder integral, employing an approach that was used in ref. [35]. 

The s>t limit of the vertical ladder is equivalent to the t>s limit of the horizontal ladder, 
and this is what we will now examine. The horizontal ladder diagram corresponds to the integral 

/L+4 \ L + 4 

Hf-T^r] (P25,mPl5,mP35,m) 1 (Pj-1 ,jPl,j,mPz,j,m) 1 {Pi,L+i,m) * (C.l) 
i=5 ^ ' i=6 

where all the internal propagators are massless while the propagators on the periphery of the 
diagram have mass m. For each of the propagators, we use (l/P) = i J °° da exp(-iaP) to 
obtain 

L+4 



ha = l' L+1 



roc ^ ^ * / 74 \ 

s L t / FT da k TT d# d lx \ TT [-4 x (C.2) 

J ° k=0 1=1 J 3=5 V m ' 

L-l L 

,l+4,m 

k=l 1=1 

Recalling that Pi 3 , m = x?- + m 2 , we express the exponent as 

-% (x T • A • x - 2 B T • x + C + m 2 a) (C.3) 
where x = (#5, • • • , xl+a)- We now integrate over x to obtain 

r°° ' - 

I La = i L+1 s L t / JJdakYldPi d-fiA- 2 exp[-i(C -B T ■ A' 1 -B + m 2 a)] . (C.4) 
J° fe=0 1=1 
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Setting x\ 2 = x\ z = x\ A = x\ x = 0, one has 



C-B T -A- 1 -B = s ^ + t^, 

A A 

Putting everything together, we have 



where s = x\ z and t = x\ 4 



(C.5) 



ha = ^ L+l 



s L t }Jda k Y[dPi d 7i A~ 2 exp 



fc=0 1=1 



..A . T 
—it— — iJ 
A 



j = s ^. + m 2 a (C.6) 



where A, D s , D t , and a are polynomials of degree L, L + 1, L + 1, and 1 in the a, (3 and 7's. 
These polynomials may be constructed using graphical rules [35,73]. 

To explore the t >> s limit of the horizontal ladder diagram, it is useful to perform the Mellin 
transform [35] 



M(e) 



dr 



' La 



with respect to r = t/m 2 . This gives 

M(e) = r(is) L (im 2 yT(l - e) F(e) 

where 



fc=0 i=l 



•iJ 



(C.7) 



(C.S 



(C.9) 



The large r behavior of i£, a is determined by the behavior of M(e) near e = 0. However, .F(e) 
diverges at e = 0, since A = n^o 01 *' so we nee d to do a Laurent expansion about e = 0. First 
we integrate by parts with respect to each of the a, to obtain 



I r°° L L ( L \ e f)L+l f p-iJ \ 

n<-^)n<ft*»(n<*j 

where we have taken e > so that the surface terms at a k = vanish. Then writing 



(CIO) 



\k=0 

we obtain 
F(e) = 



i=0 



+ 



L 1 



i=0 



+ G(e 3 ) (C.ll) 



^K,-(L + 1)K 



i=0 



L 1 



i=0 



+ 0(e 3 ) 



(C.12) 
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(a) 



(b) 



Figure 8: Factorization of the LL order contribution to the Regge limit t S> s of the horizontal L-loop 
ladder integral (a) into (2 + 2e)-dimensional bubble integrals (b). 



where we have performed (most of) the integrations over a. 



K 



'\\<l >, 



1=1 

L 



-iJ 



A l+e 

Y\ dli dai e 



a k =0 



-iJ 



K, 



poo L 

/ IT ^ ^ f2 a j~ 1 

J ° 1=1 



a k =0,k^i 
-iJ 



A l+e 



(C.13) 



a k =0,k^i,j 



The integrals in this equation may be interpreted as arising from (various deformations of) a 
horizontal L-loop ladder diagram in d = 2 + 2e dimensions [73]. The rungs with = are 
contracted to a point (and the corresponding propagator omitted), and the rungs with factors 
ea e ~ l correspond to propagators raised to the power e. Because most of the are set to zero, 
the L-loop diagram separates into a product of smaller diagrams, and the integrals factorize. For 
example, when all = 0, we have 



D„ 



i=i 



so K separates into a product of one-loop bubble diagrams (cf. fig. 8) 

K=[(i 8 y- 1 B 1 ] L 



(C.14) 



(C.15) 



where the one-loop bubble is defined as 

J i7r -ri5,m-f35,m Jo 



df3d r y 



(/3 + 7) 



l+e 



exp 



-2S 



g7 
/? + 7 



im 2 ((3 + 7) 



(C.16) 

If we are only interested in the leading log (LL) behavior as t 3> s, we may set e — >■ 0, so that 
only the first term in eq. (C.12) contributes, and B\ becomes a two-dimensional diagram. This 
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Bo 



To 



Figure 9: Integrals arising in the computation of the Regge limit of the horizontal L-loop ladder integral 
to NNLL order. The integration is in 2 + 2e dimensions. 

was noted on p. 134 of ref. [35]. (See also chapter 8 of ref. [75].) Moreover, in this limit, the 
masses of the rungs of the original ladder diagram do not affect the result. Hence, in the LL 
limit, there is no distinction between the Higgs regulator considered here and the cutoff regulator 
considered in ref. [41]. 

Taking into account the factorization of the integrals in eq. (C.12), we obtain 



F(e) = 



cL+l 



|l + [2ti + (L— l)6a - (L + 1)] + [2t 2 + t x 2 + 2(L - 2)ti&2 



with 



+(L - 2)63 + 1(L - 2)(L - 3)b 2 2 - 2Lt x - L(L - l)b 2 + X -L{L + 1)] } (C.17) 

T(l + e)\ B 2 



r(l + e) \ 7\ 
(is)* J Si' 



(isY 



to 



r(i + 6) V r 2 / r(i + e) \ 2 s 3 

(<*)« ; s?' 3 ^ («)« ; Bf 



(c.18) 



where 7\, £? 2 , T 2 , and -B 3 are the (2 + 2e)-dimensional diagrams shown in fig. 9, and defined by 
d 2+2 *x 5 1 



T 2 
£3 



= s 



i7T 1+e -Pl5 lJT t-P35,m(-P25,m) <: 



i?r 1+e 7 i?r 1 + e Pis.mPie.mPas.m^m^e) 6 
, , rf 2+2e x 5 /■ d 2+2 *x 6 1 



rd 2 ^, r<p_ 

J iir 1+t J 11 



m 1+t J m 1+e Pi^ m Pi&, m P^,mPm,m{P2^mP^y 
, , d 2+2 "x 5 f d 2+2 "x 6 f d 2+2 "x 7 I 



(C.19) 



r d 2 + 2 *x 5 r d 2+2e x 6 rd 2 

J iir 1+e J m l+t J i 



ill 1+e Pl<o t mPl&,mPl7 ,mPz5,mPm,mP37 ,m{PmP&7Y 



These integrals can be evaluated using MB techniques. 
First consider the more general integral 



/ 



d d x 5 
m d l 2 



13 



*(P 



15, m ) 



-0.1 



(P2 5 , m )- a2 (P 3 5, m )~ a3 = J gi(a u a 2 , a 3 ; d; z) u z ~ a ^ (C.20) 
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where a = a\ + a 2 + a 3 and u = m r/s and 



f , x r(-z)r(g-f + z)r(q 1 + *)r(a 3 + *) 

gia 1 ,a 2 ,a 3 ;(l;z) = — . (C.21) 

1 (ai)l (a 3 )l (2z + a) 



We will also need 

Jd 



/d cc 
— j/^iPlsT 2 (-pL5,m) ^ (P35,m) ° 3 (P35) ° 6 

= J dZ \^ly Z 92(au a 3 , a 6 ; d; (Pi,)— 23 (Pie,™) 21 (iVO* 3 «^ a+ ^ (C.22) 

where a = ai + a 3 + a 6 and 2 = z 1 + z 2 + z 3 and 

r( — z iW(— z 2 )r(— z 3 ) d 
g^a^d.Zi) = r(ai)r(fl3)r(a6) r(a - - + z) x (C.23) 

r(ai + + z 2 )T(a 3 + z 2 + z 3 )T(a 6 + z 1 + z 3 )T(d - a - a 6 - z 1 - z 3 ) 
T(a 1 + a 3 + z + z 2 )Y(d - a) 

Then we may write 



Bi = J 7WZ7^ 9i(h 0,1; 2 + 2e;^) M - 1+M 



(2tt?) 

-1-2 



Tl = / 7^^i(l,e,l;2 + 2e;z)« 

P 2 = J j^g 2 (lA,e ] 2 + 2e ] z 1 ^)g 1 (l- z u 0^- z 3] 2 + 2e-,z 4 )u- 2+ ^- z \ 

T 2 = [ 7 ^g 2 (l,l,e;2 + 2e,z l ^ 3 )g 1 (l-z 1 ,e,l-z 3 ;2 + 2e;z 4 )u- 2 - Z2 - z \ 
f d 7 z 

B? > = I m .x 7 g2(l,l,e;2 + 2e;zi ; 2, 3 )g2(l - Zi, 1 - z 3 , e; 2 + 2e; z 4 ,5,e) 

x ^(1 -2 4 ,0, 1 -z 6 ;2 + 2e;^ 7 )M~ 3+e ^ 2 ^ 5 ^ 7 . (C.24) 

Note that some of the MB representations above require an analytic continuation to e m 0. In 
some cases this may reduce the number of MB integrals. Expanding around e = (after having 
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taken the analytic continuation where required) and then expanding in small u we find 27 

+ . . . 





—2 log u + e 




log 2 u - 


r _ 


+ 6 2 


e -nv Ti = 


—2 log if + e 


"1 


log 2 u — 


—71 


+ 6 2 


2 




3 





- J log 3 u - log u + 4( 3 
3 o 



7T 2 logW + 8(3 



+ . . . 



e~ 2e ^B 2 = 41og 2 w + e 
~5 



10 i 3 4 2, 

— \0g i U+ -7T 2 logM + 4C 3 



- log 4 U + 7T 2 log 2 M + 4( 3 log U + — 7T 4 

3 to 



-2, 7ET2 = 41og 2 M + e 



-3e 7E d 



-8 log 3 u + e 



log 3 M + 27T 2 logM + 4C 3 



28 



+ . . . 

5 4 31 4 
log W H 7T 

24 & 90 



(C.25) 



+ 



— log 4 u — Atx 2 log 2 u — 16£ 3 log w 



+ 



94 38 34 

— log 5 u vr 2 log 3 u - 20( 3 log 2 u - — 7T 4 log u - 44.705 

15 9 45 



+ . . . 



Combining eqs. (C.8), (C.17), (C.18), and (C.25), we obtain 
<~2) L 



M(e) = t- 



1 + e 



^ log L+1 u + ^tt 2 log*" 1 « + (L - l)Cs log L - 2 « 



+ e 
+ 
+ 



10L 2 - 14L + 3 (L + l) 2 , L(L-2) 



180 



log^ u + 



18 
L-2 
6~ 



7l Z log^ U + 



-Calog^u 



(L + 3)(5L + 2) 7r4 bgL _ 2 u + ^ + x g2 ^ logi _ 3 u 
3o0 

(L " 2 l (L " 3) C 3 2 log L - 4 ^ + ^ 3 ) 



(C.26) 



To take the inverse Mellin transform, we use [35] 

1 



M(e) 



I = 



1 log n r 



n\ t 



(C.27) 



27 The coefficients of the n 4 terms were obtained numerically, and replaced by their probable rational equivalents. 
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so that the first three leading log terms of the ladder diagram are given by 




2 L 

lim lim I La (u, v) = —\og L v\og L u 



u + (L + 2)C 2 log 1 '" 1 u + (L — 1)C 3 hg L ~ 2 u 



>g L+2 u+ { ±±^^\og L u 

io 

^Hog L ~ 2 u + I ^{L + 1.826)7r 2 C 3 log L - 3 M 



(C.28) 



where we recall that v = m 2 /t = 1/r. This result is used in sec. 4. 
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